Improving teleportation fidelity in structured reservoirs 
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Seeking flexible methods to control quantum teleportation in open systems is an important task 
of quantum communication. In this paper, we study how the super-Ohmic, Ohmic and sub-Ohmic 
reservoirs affect teleportation of a general one-qubit state. The results revealed that the structures of 
the reservoirs play a decisive role on quality of teleportation. Particularly, the fidelity of teleportation 
may be improved by the strong backaction of the non-Markovian memory effects of the reservoir. 
The physical mechanism responsible for this improvement are determined. 
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I. INTRODUCTION 

Since the seminal work of Bennett et al. quantum 
teleportation undergoes rapid development both in the 
theory 0-0| and experiments [§-[Tl1|. Particularly, quan- 
tum teleportation via free-space channels over 97 kilome- 
tres and 143 kilometres [l3| have been experimen- 
tally realized recently by two research groups leaded by 
Pan and Zollinger, respectively. These finding rekindled 
the dream of people on constructing ultra-long-distance 
quantum communication networks on a global scale. To 
make this dream come true, one needs to face many ac- 
tual challenges. One such challenge is the decoherence of 
quantum states due to its inevitable interaction with the 
surroundings EMI. This unwanted interaction in most 
cases causes degradation of quantum correlations [l7l - [T9| , 
and therefore limits the fidelity as well as the distance of 
quantum teleportation. Even for the simplest protocol of 
one-qubit teleportation |l| , there are many processes that 
decoherence may be set in, e.g., the preparation and dis- 
tribution processes of the shared quantum channels, the 
imperfect measurements of the information sender (Al- 
ice) and receiver (Bob). All these make it significant and 
vital to investigate resistance of a teleportation protocol 
under the influence of decoherence, and this subject has 
been studied by many authors from the perspective of 
both noise channels [2fll - |26j and noisy operations [13, HH ■ 
Particularly, it has been shown that although entangle- 
ment is the necessary resource for quantum teleportation, 
in general its magnitude is not proportional to the value 
of the teleportation fidelity [HHHj]. 

In the present work, we study quantum teleportation 
of the single-qubit state in the bosonic structured reser- 
voirs of the super-Ohmic, Ohmic and sub-Ohmic types 
[29] ]. Through simulating numerically the dynamical be- 
haviors of the average teleportation fidelity, we show that 
the explicit form of the spectral density of the reservoir 
as well as the coupling strength between the channel and 
the reservoir have great impact on efficiency of quantum 
teleportation, and the quality of teleportation can be im- 



proved via the technique of reservoir engineering. 

The structure of this paper is arranged as follows. In 
Section II we recall some basic formalism related to the 
fidelity of quantum teleportation, and detailed the meth- 
ods for solving the equation of motion for the channel 
state. Then in Section III, we investigate resistance of 
quantum teleportation against decoherence induced by 
the super-Ohmic, Ohmic and sub-Ohmic reservoirs. Fi- 
nally, Section IV is devoted to a summary. 



II. THE FORMALISM 

In this work we concentrate on efficiency of quan- 
tum teleportation for a general one-qubit state, with a 
two-qubit state p serving as the quantum channel. Wc 
suppose during the teleportation process, Alice performs 
only the Bell-basis measurement, while Bob is equipped 
to perform any unitary transformation, then the maximal 
average fidelity achievable can be evaluated as [3(| 



^av = ~ + \N{p), 



(1) 



where N(p) = TrVTW, with T being the 3 x 3 positive 
matrix with elements t nm related to the Bloch sphere 
representation of p below 
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with (To being the 2x2 identity operator, 01,2,3 are the 
usual Pauli matrices, and t nm — Tr(pcr„ <g) cr TO ). 

To evaluate resistance of the teleportation protocol, we 
suppose the channel state consists of two identical qubits 
which coupled independently to their own reserv oir, with 
the single "qubit+reservoir" Hamiltonian reads [3ljj 

H = luo(t + <t- + ^ u kb\bk + ^(5fckfcC + + H.c.), (3) 

k k 

with ujo being the Bohr frequency of the two qubits, and 
the index k labels the reservoir field mode with frequency 
Wfe and the system-reservoir coupling strength . More- 
over, er_|_ (d-) are the usual Pauli raising (lowering) oper- 
ator, while b k (bk) are the bosonic creation (annihilation) 
operator. 
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Under the assumption of zero-temperature reservoir 
and nullity of initial correlation between the qubit and 
the reservoir, the reduced density matrix p s (t) for qubit 
S (S = A, B) can be obtained as [3l[ 

/pfi(0)|p(*)l a pfo(0)p(*) \ 
p S {t) = , (4) 

V p s 01 (o)p*(t) l-pfMW)] 2 J 

in the standard basis {|1), |0)}, and pfj(0) = (i\p s (0)\j). 
The density matrix p (t) depends solely on the function 
p(t), whose explicit time dependence contains the infor- 
mation on the reservoir spectral density and the coupling 
constants. To be explicitly, we consider here the family 
of reservoir spectral densities given by [29j 



J H = ^2\9k\ 2 ${u-uJk), 



(5) 



(6) 



which simplifies to 

J(w) 

in the continuum limit of reservoir modes J^k \dk\ 2 ~> 
J d(u))J(u), with uj c being the cutoff frequency and r\ the 
dimensionless coupling constant. They are related to the 
reservoir correlation time tb and the relaxation time r# 



by t b 



-1 and tr 



V 



By tuning the value of s, 



the spectral densities may be in the regime of sub-Ohmic 
(0 < s < 1), Ohmic (s = 1) or super-Ohmic (s > 1). 

For this family of reservoir spectral densities, the func- 
tion p(t) obeys the following integro-differential equation 

m 

P (t) + iu Q p{t) + I p(h)f(t - h)dh = 0, (7) 
Jo 

where the kernel function f(t—t{) = J dusJ(w)e~ luJ ( t ~ tl \ 
After a straightforward algebra, Eq. ([7]) turns into 

p(t)=p(0)-jf (two + f(t - h)dtj p{tx)dti, (8) 

In this work we take s = 1/2, 1 and 3 as three examples 
of the sub-Ohmic, Ohmic, and super-Ohmic spectral den- 
sities. The kernel function can then be integrated as @ 

( S e z), 

(« = l/2), 



[1 + iuj c {t - h)}^ 1 



tan 1 [uj c (t — 



2[l + Wc 2 (t-t 1 ) 2 ]3/4 

where si is the factorial of s, and vj = ■ 
With the help of Eq. ©, the function p(t) in Eq. © can 
be solved numerically and the two-qubit density matrix 
p{t) can be determined by the procedure of Ref. [33| . 

In fact, the equation of motion of qubit S can also be 
written in the following master equation form (by differ- 
entiating Eq. (Q} with respect to time) of the form 



P 3 (t) 



-i^[ a+ a.,p s (t)} + ^-[2a^ s (t)a + 




FIG. 1: (Color online) The \p(t)\ dependence of the average 
fidelity _F av with different mixing parameters of p(0). 



with the time-dependent decay rate T(t) and Lamb shift 
ft(t) given by [3l| 



r(t) = -2Rc 



p(t) 



p(t) 



, n(t) = -2Im 



p(t) 



(ii) 



In general, the absolute values of T(t) may be larger than 
unity. In the following, we define the normalized decay 
rate *y(t) = T(t)/T max (t) for better visualizing the corre- 
sponding phenomena, where the maximum is taken over 
the entire time region for fixed 77 and u> c . 



III. AVERAGE FIDELITY DYNAMICS 

Here we suppose the channel state is prepared initially 
in the Werner-like state I'A l| 



p(0) = r|¥><#| 



1 



(12) 



where |W) = (|00> + \ll))/V2. p(0) is entangled if 1/3 < 
r ^ 1 and separable if $5 r $C 1/3. For this kind of p(0), 
we have 



N(p) = (r + l)\p(t)\ 4 + 2(r-l)\p(t)f 



1, 



(13) 



a+a-p b (t)-p b (t)a + a-], 



(10) 



therefore the average teleportation fidelity F av is deter- 
mined by the absolute values of pit) , and this applies to 
all kinds of reservoir spectral densities with the Hamilto- 
nian model of Eq. 

By combining Eqs. (fTJ and (fT3|) , one can identify three 
different types of behavior of the \p\ dependence of F av 
(see, Fig. [1]), which depend on the values of the mixing 
parameter r: (i) for r = 1, -F av decays monotonically with 
the decrease of \p\ and arrives at the classical limiting 
value of 2/3 when \p\ = 0; (ii) for 1/3 < r < 1, F av first 
decays with the decrease of |p|, and then increase with the 
decrease of \p \ after the critical point \p\% = (1— r)/(l+r), 
and F av is larger than 2/3 only when \p\ 2 > 2|p| 2 ; (iii) for 
^ r ^ 1/3, although F av still decays with the decrease 
of \p\, and increase with the decrease of |p| after |p| > \p\ c , 
its value cannot exceed 2/3. 
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FIG. 2: (Color online) Average fidelity F av versus Dot in the 
super-Ohmic reservoir with ui c = coo, r = 1, and rj = 0.15 
(black), 0.35 (red), 0.9 (blue). The inset shows behaviors of 
the corresponding normalized decay rates 7(t). 
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FIG. 3: (Color online) Average fidelity F av versus uj^t in the 
Ohmic reservoir with uj c = ujq, r = 1, and rj = 0.3 (black), 
0.9 (red), 2.7 (blue). The inset shows behaviors of the corre- 
sponding normalized decay rates 7(t). 



In Fig. [3] we display dynamical behaviors of the aver- 
age teleportation fidelity F av for the case of super-Ohmic 
reservoir with the cutoff frequency uj c — ujq, mixing pa- 
rameter r — 1 and different coupling strengths rj. One 
can see that in the weak-coupling regime (e.g., rj = 0.15), 
-F av decays monotonically and approaches to its limiting 
value 2/3 only when uj^t — » 00. If one enlarges the cou- 
pling strength moderately (e.g., rj — 0.35), the strong 
non-Markovian effect dominates and this leads to an ob- 
vious oscillation of F av , but in the long-time limit it still 
reaches the classical limiting value 2/3. If one further 
increases the coupling strength to a much larger value 
(e.g., rj = 0.9), however, F av firstly experiences some os- 
cillations, and then approaches to a finite value larger 
than 2/3 in the long-time limit. 

To understand the physical mechanism responsible for 
the above phenomena, we plot in the inset of Fig. [5] the 
behaviors of the normalized decay rate "f(t). Clearly, for 
•q = 0.15 the decay rate maintains a finite value after one 
oscillation, and thus deprives the quantum advantage of 
the teleportation protocol in the long-time limit. For 
rj = 0.9, however, j(t) approaches zero after some oscil- 
lations due to the energy and/or information exchanging 
back and forth between the qubits and their independent 
memory reservoirs, and this leads to F av > 2/3 even in 
the long-time limit. The origin of F av (t —> 00) > 2/3 for 
rj = 0.9 can also be explained from the formation of the 
bound state in the system. In fact, for the Hamiltonian 
model of Eq. © with the spectral densities of the form 
of Eq. ([6]) , the conditions for formation of a bound state 
can be determined as 



■q > 



(s-l)lw. 

UJ 



{s e Z), 



is = 1/2), 



(14) 



by the same methodology of Ref. [32j. For the system 
parameters chosen in Fig. [5] (i.e., s = 3 and uj c = luq), 
rj = 0.9 is clearly larger than its critical value 0.5, and 




FIG. 4: (Color online) Average fidelity F& v versus uot in the 
sub-Ohmic reservoir with lo c = ujq, r = 1, and rj = 0.3 (black), 
0.55 (red), 2.1 (blue). The inset shows behaviors of the cor- 
responding normalized decay rates ~y(t). 



therefore F av (i -> 00) > 2/3. For rj = 0.15 and 0.35, 
which are smaller than 0.5 and thus F av (t —> 00) = 2/3, 
although there may be weak enhancement of F av (e.g., 
rj = 0.35) during certain uj^t regions caused by the non- 
Markovian effects. 

When considering the Ohmic (Fig. [3]) and sub-Ohmic 
(Fig. H]) reservoirs, similar phenomena are observed (the 
oscillating behaviors of F av for the Ohmic reservoir with 
rj = 0.9 and sub-Ohmic reservoir with rj — 0.55 are very 
weak and unobvious in the plots) , with however different 
critical coupling strengths for the formation of a bound 
state. This corroborates the observation derived for the 
case of super-Ohmic reservoir. 

To further understand the relations between dynami- 
cal behaviors of F av and j(t), we reexamine Eq. (|lll) . 
from which one can derive T(t) = —q(t)/q(t), with 
q(t) = \p(t)\ 2 . Therefore j(t) > when \p(t)\ decreases 
with time, and j(t) < when \p(t)\ increases with time. 
As is a monotonic increasing function of |p(i)| 2 in the 
region of \p\ 2 > \p\l = (1 — r)/(l + r), F av is increased 
whenever j(t) < 0, and decreased whenever 7(4) > 
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FIG. 5: (Color online) Comparison of the behaviors of F av 
with that of the normalized decay rate *y(t) (black) and the 
population parameter |p(t)| 2 (red) in super-Ohmic reservoir 
with u) c = uio and r\ — 0.9. The horizontal dashed line repre- 
sents constant 1/3. 



during this region. In the region of \p\ 2 < \p\ 2 , however, 
F av is a monotonic decreasing function of |p(i)| 2 , thus 
F av is decreased whenever 7(t) < 0, and increased when- 
ever j(t) > 0. But it should be note that for this case 
(i.e., \p\ 2 < \p\ 2 ), -Fav cannot be larger than its classical 
limiting value 2/3. 

In Fig. [5] we gave an exemplified plot of the time de- 
pendence of F av and "f(t) certifying explicitly the above 
general arguments. When r = 1, as mentioned above, F av 
behaves as a monotonic increasing function of |p(i)| 2 in 
the whole time region, thus it is increased when j(t) < 
and decreased when j(t) > 0. For r = 0.5 we have \p\ 2 < 
\p\l = 1/3 during the time region w i € (0.49,0.84), and 
from the solid red line shown in Fig. [5J one can see ob- 
viously that F av is decreased if j(t) < and increased if 
7(i) > 0. As the negativity of the decay rate is a signa- 
ture of strong backaction of the non-Markovian memory 



effect of the reservoir, the decrement of F av with 7(i) < 
indicates that the back flow of information from the envi- 
ronment to the channel does not always improve quality 
of teleportation. 



IV. SUMMARY 

In summary, we have studied quantum teleportation of 
a general one-qubit state in the bosonic structured reser- 
voirs with spectral densities of the super-Ohmic, Ohmic 
and sub-Ohmic types. Through comparing dynamical 
behaviors of the average teleportation fidelity F av with 
different system-reservoir parameters, we show that the 
structures of the reservoirs play a decisive role on quality 
of teleportation, and F av can be improved by the tech- 
nique of reservoir engineering, e.g., by tuning the reser- 
voir structures so that there are bound state formulated. 
We have also compared behaviors of F av with that of 
the decay rate j(t), and revealed the relationships be- 
tween the negativity/positivity of "/(t) and the enhance- 
ment/decrement of F av . 

As it is now possible to simulate and control the non- 
Markovian effect we hope our results may be cer- 
tified with state-of-the-art experiments, and at the same 
time shed some new light for exploring the relationship 
between non-Markovianity of the reservoir and the effi- 
ciency for performing certain quantum protocols. 
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